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1 ,  Introduction 

The  trustees  of  multi-employer,  jointly  trusteed  health  and  wel- 
fare funds  in  the  building  trades  face  particularly  difficult  financial 
problems  because  of  the  characteristics  of  the  industry.   The  move- 
ment of  both  employer  and  employee  from  job  to  job  and  area  to  area, 
and  the  casual  and  variable  nature  of  employment,  mean  that  a  man's 
eligibility  for  fund  benefits  cannot  be  based  solely  on  a  partic- 
ular job.  As  a  result,  a  system  of  rules  has  been  established  under 
which  a  man  becomes  eligible  for  benefits  if  he  works  a  minimum 
number  of  hours  during  a  given  period.   He  may,  by  working  the 
minimum  number  in  a  current  period,  become  eligible  for  benefits 
at  a  future  date,  or  he  may  have  earned  current  benefits  because 
of  hours  he  has  worked  in  the  past.  The  minimum  number  of  hours, 
the  period  during  which  they  must  be  worked,  and  the  length  of 
time  a  worker  remains  eligible,  vary  greatly  from  fund  to  fund. 

Essentially,  ti.e  ;;oals  of  fund  trustees  are  to  provide 
benefits  for  the  "workers  active  at  the  trade,"  to  use  the  maxi- 
mum Sums  available  for  current  benefits,  and  at  the  same  time  to 
maintain  fund  solvency.   Their  decisions  concerning  eligibility 
rules,  the  amount  of  money  that  will  be  spent  for  benefits  per 
eligible  worker,  and  the  sum  of  money  to  be  held  in  reserve 
determine  the  fund's  financial  position.  Because  of  the  fluc- 
tuations in  employment  in  the  building  trades,  however,  the 
trustees  can  never  predict  in  advance  just  how  many  men  will 
be  working,  the  number  of  hours  they  will  work,  and  how 
many  of  them  will  be  achieving  eligibility.   Consequently, 
they  can  never  be  sure  exactly  what  their  income  or  expenses 
will  be  at  any  one  time. 

The  fund's  financial  problems  atfect  the  degree  to  which 
the  trustees'  will  achieve  their  goals.   When  increased  expense, 
reduced  income,  or  inadequate  reserves  threaten  the  fund  with 
incolvency,  the  trustees  can  retrieve  the  fund's  position  by 
changing  the  eligibility  rules  to  reduce  the  number  of  benefi- 
ciaries, by  shortening  the  period  during  which  members  remain 
eligible,  or  by  cutting  down  on  the  benefits  available.  All  of 
these  solutions  have  been  used  by  various  funds,  but  each  one 
reduces  the  effectiveness  of  the  health  and  welfare  plan. 
Changes  in  eligibility  rules  and  reduction  in  benefits  under- 
mine the  employees'  confidence  in  the  fund,  and  may  lead  to 
their  purchase  of  coverage  outside  the  fund.   This  in  turn  may 
cause  unrest  in  the  union  and  internal  political  difficulties. 
Even  when  a  plan  is  running  successfully,  many  members  prefer 
to  receive  the  amounts  contributed  by  enployers  in  cash  rather 
than  benefits.  When  their  confidence  in  the  benefit  plan  is 
undermined  by  a  series  of  changes  in  eligibility  and  benefits, 
they  will  certainly  prefer  the  cash.   Frequent  changes  may, 
in  fact,  lead  to  the  termination  of  the  plan. 


In  rather  simplified  terms,  the  essence  of  fund  management  lies 
in  control  of  the  money  flowing  into  and  out  of  the  fund.   The  flow 
of  income  into  the  fund  during  a  particular  period  is  deteruined 
largely  by  the  contribution  rate  and  the  working  hours  to  which  the 
rate  applies.   It  is  also,  of  course^  affected  by  earnings  from 
investment  of  reserves  and  by  reciprocal  agreements  between  funds. 

Money  flov;s  out  of  the  fund  primarily  for  the  purchase  or  pro- 
vision of  benefits  and   for  operating  expenses.   The  amounts  spent  on 
benefits,  assuming  tiiat  tb.ese  are  purchased  rather  than  provided 
directly,  depend  on  the  premium  rate  and  the  number  of  eligible 
workers.   This  last,  in  turn,  depends  upon  the  number  of  members 
who  work  the  required  number  of  hours  called  for  by  the  eligibility 
rules.   The  operating  expenses  of  the  fund  are  determined  by  the 
size  of  the  fund,  the  type  of  administration,  the  efficiency  of 
the  administrator,  and  the  range  of  functions  performed. 

The  timing  of  flop's  of  money  into  and  out  of  the  fund  is  an 
important  determinai.l  of:    the  fund's  financial  position.   Since 
the  eligibility  rules  ,,c:.ieraHy  provide  for  future  benefit  cover- 
age on  the  basis  of  hourc  vrarked  in  a  current  period,  theoreti- 
cally current  income  pays  for  future  benefits.   The  timing  of 
money  outflows  is  affected  not  only  by  the  number  of  workers  who 
become  eligible,  but  by  the  length  of  time  over  which  coverage 
is  continued.   The  eligibility  rule  and  the  period  over  which 
benefits  are  payable  create  a  future  liability  for  the  trustees, 
even  though  their  legal  liability  continues  only  so  long  as 
there  is  money  in  the  fund. 

Not  infrequently,  more  money  is  flowing  out  of  the  fund  than 
is  being  replaced  by  income  into  it.   For  this  reason  a  reserve 
sum  of  money  is  needed  to  serve  as  a  reservoir  which  can  be 
tapped  to  meet  expenses,  and  to  maintain  solvency.   The  reservoir 
may  also  enable  the  trustees  to  adjust  eligibility  rules  during 
periods  of  depressed  employment  conditions  to  provide  for  the 
continued  coverage  of  workers  active  in  the  trade  but  not 
currently  employed.   In  some  funds  this  reservoir  is  created 
by  provisions  specifying  that  benefits  will  begin  only  after 
the  fund  has  been  collecting  income  for  some  period  of  time. 

Inevitably,  fund  trustees  find  themselves  in  a  sort  of 
three-cornered  tug-of-war  as  they  seek  a  mi.xture  of  eligibility 
rules  and  benefit  package  which  will  balance  their  conflicting 
goals: 

1.  To  provide  the  maximum  amount  of  benefits  possible. per 
eligible  employee; 

2.  To  cover  as  many  active  workers  in  the  trade  as  possible; 

3.  To  ensure  that  the  fund  remains  solvent. 


The  sources  of  conflict  among  the  goals  --wc-   obvious.   It  the  eligi- 
bility rule  is  very  stringent^  few  employees  x^ill  be  covered  but 
the  benefits  per  covered  employee  can  be  gRncvous,  Alternatively^ 
if  the  eligibility  rule  is  not  made  more  restrictive^  and  yet  the 
benefits  offered  eligible  employees  are  expanded^  the  fund  runs  a 
higher  risk  of  insolvency. 

If  the  trustees  are  to  trade  off  among  these  three  goals  most 
effectively,  they  must  have  a  clear  understanding  of  the  way  the 
elements  involved  interact.   If,  for  exam.ple_,  workers  are  required 
to  work  ten  more  hours  to  become  eligible  for  benefits^  exactly  how 
many  employees  would  then  fail  to  achieve  eligibility?  And  how 
much  would  the  odds  of  the  fund's  becoming  insolvent  be  decreased? 
In  sum,  what  the  trustees  need  is  a  clear  idea  of  the  quantitative 
impact  of  a  variety  of  combinations  of  eligibility  rules  and 
benefit  costs.   They  need  to  know  not  only  whether  or  not  a 
change  in  expenditures  for  benefits  or  in  the  eligibility  rule 
will  achieve  one  of  their  goals  to  a  greater  or  lesser  degree,  but 
they  need  also  to  know  how  much  of  each  of  their  goals  would  be 
achieved  or  sacrificed  were  they  to  make  the  change.  As  matters 
now  stand,  however,  this  sort  of  careful  analysis  of  the  inter- 
actions is  not  ordinarily  undertaken,  and  as  a  result  many  health 
and  welfare  plans  are  not  meeting  as  well  as  they  could  and 
should  the  goals  for  which  they  v-'c-rfi  established. 

Even  though  the  trus'.c-cf  c-.i'i  .^•ik','  r,   multiplicity  of  choices 
among  combinations  of  eligibility  rules,  expenditures  for  benefits, 
and  fund  reserve  requirements,  some  aspects  of  the  way  these 
interact  are  evident.   Increases  in  benefit  expense,  for  example, 
will  always  mean  increased  expenditures  for  the  fund  if  all  olli'^r 
aspects  of  fund  operation  remain  the  same.  And,  given  the  usu:;! 
spread  of  working  hours  among  members,  a  reduction  in  the  aunibrr 
of  hours  needed  for  eligibility  will  alf-o  ^vcr,n   increased  expendi- 
tures because  there  will  be  a  larger  numbpi  o!  eligible  workers. 
Because  of  the  great  number  of  possible  combiny t ions ,  and  the 
wide  variety  of  possible  results,  it  is  almo-:-f  i '"■v.-'islb  l,c  £or 
the  trustees  to  know  ex??ctly  how  one  change,  oi  several  changes, 
will  affect  all  tl\e  o'.hor,  o'pccls  of  fund  operation. 

The  simulation  model  described  in  this  paper  attempts  to 
depict  the  essentials  of  fund  behavior  in  such  a  fashion  that 
the  interactions  described  above  may  be  determined.  In  order 
to  carry  out  the  analysis  the  following  information  is  required i 

1,  The  number  of  m.en  active  at  the  trade  in  the  jurisdic- 
tionj 

2.  A  history  ever  time  of  the  number  of  hours  worked  by 
the  members s 
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3c  The  eligibility  rule^ 

4„  The  contribution  rate| 

5o  The  premium  ratej 

6.  Administrative  expenses^ 

7.  The  aiTxaunt  of  the  fund  reserve  at  the  time  of  analysis. 


1.1  Scope  of  Analysis 

Of  the  ma,ny  paramecers  which  deterTiir":  ch2  behavior  of  the  eco- 
nomic reserves  in  the  fund^  the  trustee-  la  I'act  control  only  the 
parameter  (X^L)  of  the  eligibility  rule^  the  initial  reserves  U  in 
the  fund  before  benefit  payments  begin,  the  insurance  premium  r  per 
member  covered,  and  to  some  extent  the  employer  contribution  rate  k. 
Our  objective,  which  we  will  restate  in  slightly  different  language 
in  1.4,  is  to  determine  the  way  the  probability  of  insolvency,  and 
the  mean  and  variance  of  members  covered  per  time  period  vary  with 
changes  in  these  controllable  parameters. 

In  Part  I  of  this  paper  we  describe  a  model  of  the  random  process 

generating  hours  worked  per  time  period  by  each  member  active  at  the 

trade,  and  then  show  how  the  behavior  of  eligibility  rules  may  be 

investigated  using  simulation  techniques.   The  class  of  eligibility 

rules  we  examine  in  some  detail  may  be  described  as  follows: 

A  member  must  work  at  least  X  hours  during 
time  periods  t-^L  through  t-1  in  order  to  be 
covered  during  period  t;  t=:l,2,,,„  and  L^  1. 

In  Part  II  we  present  a  series  of  graphs  that  describe  how  the 
probability  of  insolvency,  the  mean  number  of  members  covered  per 
period  and  the  variance  of  this  mean  vary  with  changes  in  values  of 
controllable  parameters.   These  graphs  are  derived  by  a  mixture  of 
mathematical  analysis  and  monte  carlo  simulation,  as  the  reader 
will  subsequently  see.   Since  (X,L)  is  the  only  controllable  param- 
eter influencing  the  mean  and  variance  of  members  covered  per 


period,  we  display  their  behavior  as  a  function  of  (X,L)  independent 
of  graphs  of  the  probability  of  insolvency  as  a  function  of  control- 
lable parameters; 

MEAN  OF  NUMBER 

OF  MEMBERS  COVERED  PER  ;'ERIOD 


L" 

-    L' 

L 
->  X 
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One  such  graph  is  displayed  for  each  (X,L)  pair. 


The  mathematical  analysis  underlying  these  graphs  is  presented 
in  Part  1,   We  define  symbols  in  1.2,  set  down  the  assumptions 
underlying  the  model  generating  "hours  worked  by  the  ith  worker 
in  period  t"  in  1<>3,  and  restate  objectives  in  more  formal  lan- 
guage in  lo4„   Section  Io5  defines  auxiliary  symbols  used  in  the 
following  sections. 

In  section  2«1  we  derive  the  mean  vector  and  covariance  matrix 

of  a  random  vector  &  ,  and  use  them  to  determine  the  mean  and 

— t' 

variance  in  any  given  period  of  the  amount  of  economic  reservfs  if 
the  fund.   The  core  of  the  sim'alation  routine  consists  of  gener- 
ating values  of  6  and  then  using  these  values  to  determine  how 
the  probability  of  insolvency  varies  with  changes  in  parameter 
values.  We  also  present  the  mean  and  variance  of  the  average 
number  of  members  covered  per  time  period  as  a  function  of  quan- 
tities calculated  in  2.4  and  2.5,  A  mutivariate  r<r.ntral  limit 
theorem  is  applied  to  8   in  2„2,  demonstrating  that  6   is  asymp- 
totically  Multinormalj  we  exploit  this  fact  later j,  and  by  doing 
so,  greatly  reduce  computation  time.   Sections  2.3,  2.4,  and  2.5 
show  how  to  calculate  elements  of  the  mean  vector  and  covariance 
matrix  of  6  . 

In  section  3  we  describe  the  structure  of  the  simulation 
routine.  Section  3.1  shows  how  the  periodicity  of  elements  of  the 
covariance  matrix  of  6  may  be  exploited  to  enable  us  to  find  the 
square  root  S  of  this  matrix,  no  matter  how  large  its  order  may  be. 
(The  matrix  S  is  used  to  transform,  sequence  of  standardized  normal 
random  deviates  into  values  of  6  ).   Finally,  in  section  3.2  we 


outline  the  computer  program  used  In  the  simulation. 

Part  II  applies  computer  simulation  routine  and  mathematical 

results  of  Part  I  to  data  from  the  a  building  trades  union.   Results 

are  presented  in  the  form  of  the  graphs  described  earlier. 


1.2  Definitions 


Define 


h   -  number  of  hours  worked  by  the  ith  worker  in  period  t, 

I  -  income  in  the  tth  period, 

a  -  expense  in  the  tth  period, 

k  -  employer  contribution  rate    per     period, 

n  -  number  of  employ<ej  covered  in  the  tth  period, 

N  -  number  of  members  active  at  the  trade, 

U  -  amount  of  reserves  in  the  fund  at  the  end  of  period  t, 

r  -  insurance  premium  per  covered  enf)loyee, 

L   -  a  parameter  of  the  eligibility  rule;  the  number  of  past 
periods  over  which  hours  are  accrued  in  determining 
a  worker's  eligibility  for  any  given  period, 

X   -  number  of  hours  a  worker  must  work  during  periods  t-L 
through  t- 1  to  be  covered  in  the  tth  period. 


t 
U.  =  (U^  -  ta)  +  Z  (I   -  rn  )   ,  (la) 

T  =  i 

where 

N 

r  =  2  »<  h    •  (lb) 

^   i=l   ^^ 

Also  dafine 

n  =  (a,  k,  r,  U^,  N,  X,  L)    , 

p  (t)=   1  -  F(U  >  Oj  •  '     :  ;)   •  t:b?  probability  ot  insolvency 
n  1,  ^        ■.     _ 

wrtnxn  the  tirre  interval  .'  i.G  t  giv^^n  ,,  . 

E  (t)^   tha  excected  vaiting  tir^:=  to  inficlvency  given  p^  ccndicionrii 
on  intclvency  occurring  within  the  tirre  interval  0  to  t. 


Then 


-  10 


1.3  Assumptions 

I  The  eligibility  rule  takes  the  form;  a  worker  must  work  at  least 
X  hours  during  periods  t-L  through  t-1  to  be  covered  in  period  t, 
t=l,2,...  and  L  >  1. 


II    The  behavior  of  \\.      as  t  increases  may  be  represented  as 
h.^  =  f  (t)-+  e.^   ,  all  i  ,   t=l,2,... 

where  the  random  variables  {e.  ,    i=l,2^...,N;  t=l,2,...3 

are  mutually  independent  and  identically  distributed  with 

2 
mean  0  and  variance  a   ,    and  f(t)  is  some  function  of  time, 


III    The  number  N  of  members  active  at  the  trade  is  constant  over 
time  periods  t=lj2,... 


1,4  Objectives 

Given  the  ordered  7-tuple  (a,  k^r,  U  ,  N,  X,  L)  £  n  and  Assumptions 

I,    II J   and  III  what  is: 

(a)  the  expected  waiting  time  E  (t)j 

(b)  the  probability  P^Ct); 

(c)  the  marginal  probability  distribution  of  U  {Pi  :■'.    tj 

(d)  how  do  (a),  (b),  and  (c)  vary  with  changes  in  the  elements  of  n? 

(e)  how  do  the  mean  and  variance  of  n  vary  with  changes  in  the 
elements  of  Q? 
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1.5  Further  Definitions 

It  will  be  useful  to  re-express  (la)in  the  form 

U,  =  Uq  +  Z     5^  ,  (2) 

T  =  i 

where 

*%/  '*v  'N^  »  _  . 

^-"t-Vi     '  <^^ 

Define  for  t=1^2,,,,  and  for  t=1,2j... 

E(§j.)  =  6^   >    V(S|.)  =  Cg    ,    Cov(S|.,6^)  =  0^^,  t;^  T, 

t-1  t-1 

X  =  ;:  -  f(t)  ,   p.  =  P(   E  €.   >  X  -  f(t)  )  =  P(   Z  G.   >  X  )   ,    all  i, 

T=t-L  ^^  T=t-L  ^"^ 


i;^  =  -a  +  f(t)  -  f(t-l)   , 

6^  =  (6^  6-  ...  6   ...  b)^ 

~t  12       T       t 


2.    Initial  Results 

2 . 1  Mean  Vector  and  Variance-Covariance  Matrix  of  5 


We  will  show  that  the  random  vector  6   is  distributed  with  mean 


where 


6^  =  (6  6   ...  6  )" 

— t        1    ^         T 


6  =  t  -  rNp   ,   T=.,2,...,t   ,  (4) 

T      T         T 


'v.e  v;ill  ii.oulye  in  abuse  of  flotation  and  give  the  symbol  t  dual  meaning. 
U'hen  it  appears  in  a  subscript  it  denotes  time  period  t,   As  a  super- 
script it  denotes  the  transpose  of  a  matrix  or  vector.   The  context 
will  make  the  meaning  clear. 
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and  with  symmetric  covariance  matrix 


0-.   »  •  o  o  ^       •  ♦  I 

5^      It 


•2 
0  ,  •  •  •  O-      •  o 

t1     o 

T 


0   .  »  »  0   ,  ,  , 

tl       tT 


It 


Tt 


•2 

"6. 


(5) 


where  the  elements  of  E  are  determined  as  follows:   defining  for  1  <  i  <  L 
and  i=l,2, , . . ^N, 


%,T+i  =  ^^^,T+i  ^T^ 


(6a) 


and 


T+i-1  ^ 

q    _  =  P(  E    e.   >  X    ,    E  €.   >  X  )   , 


T-1 

2 

t=T-L 


(6b) 


we  have 


Ti^x'+i 


if 


and 


TT       1   2 

'^  N[q   j_«  -  P  P  _,J  -  rkNc  T=T'+i 

■■      T,T+i  T     T+i  T,T+i 


2       2  2  2 

0   =  N{r  p  (1-p  )  +  k  0  )  ,        T=1^2y...,t, 
o        T    T      e 


(6c) 


(7) 


From  (2),  (5),  (6a),  (6b),  (6c),  and  (7)  it  follows  that 


E(U  )  =  U-  -  at  +  kf (t)  -  rN  E  P   , 

t       U  ,    T 

T  =  i 


(8) 


t  ^      t        t-1   t 
V(U^)  =V(E5)=Eac.+2E    E   c^,. 


T=l         T=l    T       T  =  l    t'>T 


TT 


(9) 


i 


See  formula  (10)  for  a  definition  of  the  random  variables  y.   ,.  • 
^      ^  ■'x,  r+£ 
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We  may  calculate  both  the  mean  and  variance  of  the  number  of  workers 
covered  per  time  period  as  functions  of  p  ,  t=1,2,,.,  and  q    ..,  1  <  i  <  L, 
and  so  facilitate  analysis  of  their  behavior  as  functions  of  the  parameter 
(X;L)  of  the  eligibility  rule:   Define 

1   t  ^ 
n     =  —  Z,  n   -  average  number  of  members  covered  per  time  period 

't     t  T  =  l    T  ^  f  f 

over  t  periods  of  time. 
We  show  below  that 


-1        T 


and 


v<\>  =  -2  ii  ^<^-^>  +^  T?i  ^  ^V' '  p/t'^  '   ^^°^^ 

t  t      t'>t 


]'i 


PROOFS ;   Before  proving  (4)  through  i^) ,   observe  that  the  data  generating 
process  which  gives  n  may  be  thought  oi  as  a  Bernoulli  process:   define 


for  given  L 


y. 


i,t+l 


'1        ^   i',.,  Il  - 
T=t-L+1 

if  all  i,   all  t. 


-  Z    h.  <  X 
T=t-L+1  ^^ 


If  we  define 


t 

E    e 


X  =  X  -  f(t)  and  s.   = 


then  we  may  write  (11)  as 


^i,t+l 


s.  >  X^ 

IT  -  t 


s.   <  X 
it    t 


,   all  i,  all  t. 


(H) 


by  virtue  of  II. 

Since  the  €,  s  are  mutually  independent  by  Assumption  11^  so  are 

the  random  variables  y.  ,  and  furthermore,  since 

•'it-'  ' 

P(s.^  >  X^)  =  P(s^^  >  Xj.)  s  p^   ,   all  1  <  i,  j  <  N;   all  t,      (12) 


for  a  given  t  we  may  regard  y,  ^  y«  ^ •  •  • ^  y  ,»..,   y       as  values  generated 

by  a  Bernoulli  process  with  parameter  p  . 

It  will  be  convenient  to  work  with  the  y  s  rather  than  n  in  the 

■^it  t 

subsequent  proofs; 
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Proof  of  (4) ;   From  (1)  and  (3)  we  have  for  t  >  1, 

t   N 

5^  =  i^n  -  cct  +   f(t)  +  E   Z  (k  e   -  ry.  )] 
t     u  1  .  ,     It     IT 

T=l  l.  =  l 

t-1  N 
-  [Un  ~  ccit-l)   +   f(t-l)  +  E   E  (k  €   -  ry  )] 

T=l  1=1 

N 
=  -a  +  f(t)  -  f(t-l)  +  E   (k  e.j.  -  ry^^)  . 

i=l 

Using  II  and  the  fact  that  ECy^^^^)  =  p^  all  i, 


6^  =  -a  +  f(t)  -  f (t-1)  -  rNp^. 


Proof  of  (7) ;  Assumption  II  implies  that  y.   is  independent  of  y ,  ^ 
j^i^  and  also  that  y.   is  independent  of  e .  ^  j^i.   Thus,  using  II  and  (12), 

N  N   -,  . 


V(6  )  =  V(  E  (ke 
■"      i  1    "■ 


T  -  '^i?    ^  %f/  '^'iT^  ^  '   '^>'iT> 


2  2  2 

=  Nr  p  (1-p  )  +  Nk  a  ,   T=l,2,.o.,t„ 

T  ■     T  6 

for  Assumptions  I  and  II  together  im.plv  that  y,   and  e.   are  independent. 

It      i  t 

Proof  of  (6) ;   To  establish  (6)  observe  that  I  and  II  imply  that 


(a)  when  t  ^  t'  +  i,  i=l,2, , « „  ,L;,   y.   is  independent  of  e  .  , 
for  all  i  and  jj 

(b)  when  t  =  t°  +  £   for  some  i,  1  <  i  <  L,   y.   is  independent  of 

'   —   —  '    It 

e.  J  only  when  il   j,  but  y.   and  e.  ,  are  correlated 5 
jT     ■'  -^^  ■'It      it  ' 

(c)  when  i^j,  y.   and  y.  ,  are  independent  for  all  t  and  t  j 

^iT 

(d)  when  i=j  and  T^T'.Aand  y.  ,  are  correlated  if  T=T'+i-l  and 

/\       JT 

1  <  i  <  L^  and  are  independent  otherwise. 
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An  example  will  help  to  clarify  the  meaning  of  (a)  ^  (b) ,  and  (d)  .   •■  ippose 

I  =  3  and  we  displcy  y    and  e.  ,  as  shown  below  for  t'=1^2,...   An  x  in  the 

i-r       3-T 

(r')'cl\   column  and  -|  th  row  indicates  that  y.   and  e.  ,  are  correlated. 

t':   I     2      3      4       5     6     7      .      . 

e 


'il     i2 


^^13 


i4 


€     e     e 
i5    i6    i7 


2  ^i.. 


^i: 


y-: 


'      ^5 


'      ^'iG 


X 


X 


X 


8   y 


i8 


X 


For  example,  when  t  ^  t'  +1,  t'  +2,  t'  +  3,  y.   and  e.  ,  are  uncorrelated . 

IT         IT 

while  if,  say,  t  =  t  +2,  y   and  e   ,  are  correlated.   Furthermore  since 

iT      It 

y.   ,,  _ ,  say,  is  correlated  with  €.   ,  ,,, ,  e.   ,  . ,  ,  and  e.   ,  ,  and  y.   ,  ,, 
i,T  +3'    ■' '  i,t'+Z'      i,t  +1^      ^f-r  ^i,T  -!^: 

is  correlated  with  e.   ,     e.  ..,    e.   .  , ,  then  y.      and  y.      are 

1,T  +3'    1,T  +2-'    l^T  +1  l^T  -Kt         1,T  -1-3 

correlated.   However,  y.   ,  ,  is  uncorrelated  with  y.   ,  ,. 

'l,T+4  1,T+1 
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Case  I;  t=t' H  ,      1  £.  ;.  1  L 


1.   First  we  show  that  we  must  evaluate  the  term  E(5  6  ,  ) 

T  T+Ji 

In  order     to  prove    (6)   when   t=t'+^   and   1  <  ^  <   L: 

Cov(6  5        )    =  E(\b      -   6    i[    6  _^     -   6        ^) 


=  E(S  &        )    -   6  6 


\;e  now  show  that  for  all  i 


Cov(y 


,  e.  )  =  E(y.     e.  ) 
where  b  is  defined  in  (15)  below  and  that 

for  we  need  these  results  to  evaluate  (14) «  We  first  evaluate  E(y.  4.0   ^j    )• 
Remember  that 

"'^^^   (_0      T'=T4.i-L  ^^   <X^^^ 
We  have  assumed  that  r+i-l  >  t'  >  T+i-L,  so  for  notational  convenience  define 

T+^-1 

b  =  (    E    e   ,)  -  €      all  i. 
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Now  given  b  =  b,  and  e.   --  z, 

'  XT  ' 

'^\-r^,   =  I  I  b,  .)  =  j   if   bfz 

r+£ 
and 

Defining  for  all  i  and  -r , 


F  (z)  ^  P(€,   <  z)  and  F,  (b)  =  P(b  <  b) 
e  in  b   ■ 

we  have 

CO     00 

-co  ""OO 


09        00 


r     r  z  dF^(z)dF^(u)  (15a) 


(x^,-)  - 


This  last  fonuula  will  be  used  to  evaluate  c  ^.   numerically. 
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We 


may  evaluate  E(y,   ..Y,  )  in  a  similar  fashion.  Define 


T-1 
Z 

t=T+^-L 


V^^  =     ^    ^it   •  ^^^^^ 


T+i-1 
t  =  T 

and 

T+^-L-1 

?i,  =     I    ^i,   ,  (15d) 

t=T-L 
so  that 

>  X 

~   T 


<%*  Arf 


y   =  •<    if  r.   +  V. 

iT^      I  IT      IT 

.0  <  X 

T 

and 


Hence  if  r,   =  P,   V,   =  V,  and  w.     =  W,  then 
iT        It    '      i,T+i 

rl  >  X 

P(y.  =  i|r,v)  =4   if  r,  +  V. 

It     '        /  IT    IT 

^-0  <  X 

T 


^(^i/T+i  =  Mr,  w)  .  j   if   v^^fw^^^^^ 


-   T+i 


T+i 


Define  for  all  i  and  t,   and  for  1  <  i  <  L 

T 

Fw   (")  =  1  -  %        (w)  =  P(W,     <  w)  , 
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and 

F^  (g)  =  1  -  G^  (g)  =  p(r^^  <  g)  , 

T  T 

Thus  using  the  fact  that  V.  ,  W.  .„    ,   and  r.   are  mutually  independent. 

since  the  e.  's  are  mutually  independent  for  all  i  and  t,  •  i  hav 

'  )  ;ad   a  formula  for  calculating  q    ,.  when  l<  I  <   L: 

T,T+i         —    — 

00 

^       T  T+i  T 

-co 

3«   We  evaluate  E(&  5   )  by  noting  that  from  (3)^ 

T  T* * 

N  N 

N  N 

+  (  Z  ke.   -  ry.  )  (  Z  kc.     -  ry.   , J 
i=l   ^"^  i.=l   ^^'^■*''^     i.^T+i'^ 

where 

^^  =  -a  +  f (t)  -  f (t-1)     , 

i^+^   =  -a  +  f(T+i)  -  f(t+i-l)  , 
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so  that  using  (12)  and  II, 


N  N 


+  E(  E  (ke.   -  ry.  )  E  (ke .     -  ry.  ^, 

i^j^     IT      'it   ^^j^     J^T+i      ^J^T+i 

To  finish  proof  of  Case  I  we  show  that 

By  virtue  of  (a),  (b) ,  i'^) ;   ^^"^   (<^)  stated  at  the  outset  of  the  proof 
this  reduces  to 

''  i=l  ^^'-'i>T+i^  -^  r^NCN-Dp^P^^^ 

N 
-rk  E  E(y     e  ). 

j^l      J,T+i  JT 

Using  the  definitions  of  q        ,„,   P     and  c        ,.    the  above  may  be  written  as 


r^Nq^      . .  +  r^N(N-l)   p  p        -   rkNc 

T>T+*  T    T+i  T,T+i 


From  (14),    (16),   and   (17),   we  have 


Cov(6   ,B    ^,.>    =   ^^r^rn   ^^^^t+Zt  "   "^^T^T+i^ 

+  ^^\,,^,   +  r'N(N-l)p^P^^^   -   rkNc^^^^p 
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T    T+i 
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Since 

^  ^  a.*   =   (5      -   Np   )a    ,,    -   Np        ) 

T    T+^  *T  T     ^T+i  T+i 


=    C,^,,.,     -     -^?,i,^,    -     -^P,^,^,    +    ^'^'PT^T+i  ' 


*\/  /\^ 


2  2 

Cov(6    ,6    , J    =   r   Nq         , „   "    r  Np  p   ^^   -   rkNc 

which  completes  the  proof  for  the  case  t  =T'+i,   1  <  i  <  L. 


Case  II:  Ti^T'+^,  1  <  -«  <  L^ 


If  r^r'+Z     we  have 


*Srf     'N^ 


Cov(6  ,6  ,)  =  E(5  6  ,)  -  6  6  , 

T   T  T  T         T  T 

as  before,  and  by  analogy  with  (16)  we  have 

E(6  6  ,)  =  ^  ^  ,  -  rNp  ^  ,  -  rNp  ,^ 

T  T         T^T  T  T  T 

N  N 

+  E([  E  (k?   -  ry  )]  ]  E  (ke   ,  -  ry   )])   . 
i=l   ^      ^      j=l   J       J^ 

2  2 
The  expectation  on  the  right  hand  side  above  reduces  to  r  N  P  p  ,  by 

T  T 

use  of  (a),  (b),  and  (c),  giving 

Cov(6^,6^.)  =  C^C^.  -  rN(p^^^  -  rNp^p^,)  -  5^6^,  , 
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but  since 

&  6  ,  =  (;  -  rNp  )  (;  ,  -  rNp  ,) 

T  T       *T         T      T  T 


=  i    i    ,    -    '^NCp  i;  ,  +  P  ,^   -  rNp  p  ,)   , 


Cov(6  ,6  ,)  =  0   ,    Ti^T'+i  ,     1  <  i  <  L. 

T   T  —    — 


Proof  of  (10);   The  proof  of  (10a)  and  (10b)  is  in  fact  imbedded  in  that  of 
(9),  for  from  the  proof  of  (9)  we  have 

v(\)  -  v(  i  ii  \>  =  ;2  ^<  4i  \> 


'>^  *%* 


1   t    ~     2   t-1   t 

-o  E,  V(n  )  +  -o  Z.   E   Cov(n  ,n  ,) 

t      ^=^    ^    t^  ^=^  t'>t      ^  "^ 


N   t     .,    .  .  2N  t;:l   t 

tq   ... 

TT       t  T 


"2  xil  Pt(^"Pt>  "'T  t=1   S   ^''-'  '  P-''-'^ 


t'>t 
Formula  (10a)  is  obvious. 
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2,2  Asymptotic  Normality  of  6 


The  following  Lemma  is  of  conpiderable  practical  importance  because 
it  allows  us  to  simulate  a  sample  realization  (S, ;6„ ^ . .  •  ;5  j  by  generating 
only  t  random  Normal  deviates  when  N  is  large,  in  place  of  generating  2Nt 
random  numbers:   N  ys  and  N  es  for  each  of  t  6s. 

Lemma  1:   As  N  -♦>  »  the  random  vector  5   is  asymptotically  distributed 
as  a  multivariate  Normal  vector  with  mean  6  and  covariance  matrix  Z. 


PROOF : 


Define  the  2t  x  1  vector 


and  the  (t  x  2t)  matrix 


(18) 


A  = 


k  0  0  0...0  r  0  0 
OkOO,..  GOrO 
OOkO...   GOOr 


G    G    0    G 


0  ...  0  G  0 
G  ...  0  0  G 
G   ...   G    G    0 


G    G    G    G 


G    G 


(19) 
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so  that 

N 

Z  A  X.  =  6    .  (20) 

If  we  consider  x.  as  a  random  vector,  observe  that  Assumption  II  inqilies 

that  the  {x.,  i=l,2,»..,N)  are  mutually  independent  and  identically  distributed 

with  a  mean  x  and  a  symmetric  covariance  matrix  T« 

* 

We  may  now  use  a  multivariate  central  limit  theorem  to  prove  that  the 

random  vector  6   is  asymptotically  multivariate  Normal  with  mean  vector 

~  t 

N  A  X  and  covariance  matrix  N  A  T  A  as  N  -►  »: 

Theorem; 

Let  the  2t  component  random  vectors  x^  t)e  independent  and  identically 

distributed  with  means  x  and  covariance  matrices  E(x.  -  x)  (x-  -  x)  =  Z» 

~    ^   N  ^    _^        ^ 

Then  the  limiting  distribution  of  Z  =  J E  (x.  -  x)  as  N  -  «  is 

~  n/  N  i=l   ^ 

£^1''\l\9.4)    =   (2n)   ^e  ^"  =   "  |T|"^   . 


The  theorem  thus  in^lies  that  as  N  -►  oo, 

Xs  Z   X.  =  VNz  +  NX  ~  f^^(X.  |Nx,NT) 
i=l   ^  N         _ 


or 


^    N         rt)  ~   ~      t 

6^  =  Z  A  X.   ~  f  „  (B^iNAx  ,  NAT  A  ) 
~t    .,    ~  ~X  N   ~t '  =r-  '  ^1^  — 


* 

See  T,  W.  Anderson,  In:roduction  to  Multivariate  Statistical  Analysis 

(John  Wiley  &  Sons,  New  York,  1958)  . 
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2.3  Structure  of  the  Covariance  Matrices  T  and  Z 


Th 


e  covariance  matrix  T  of  x. ^  all  i,  may  be  written  as 


P  I 


H 


(21a) 


where  P  =  cr   ^  the  matrix  I  is  a  (t  x  t)  identity  matrix,  and  the  (t  x  t) 
matrices  H  and  Y  are  defined  below. 

Before  we  define  H  and  Y  explicitly,  observe  that  since  the  matrix  A 
of  (19)  may  be  written  as  (k  I  ^   1) ,   where  I  is  identity  matrix  of  order  t, 
and  since  Z  =  A  T  A  ,  we  may  use  (2.1)  to  represent  Z  in  the  following 
convenient  form: 


2  2 

Z  =   pk  I   +  2rk<I>  +  r  Y 


(21b) 


where 


for 


*  £  H     +  H  , 


Z=ATA=(kI        rl) 


pi  H 


H 


kl 


rl 


=  pk  I  +  rk(H*^  +  H)    +  r^y 

Both  H  and  Y  are  determined  by  the  structure  of  the  eligibility  rule. 

Given  an  eligibility  rule  as  stated  in  Assumption  I,  a  worker  must  work  a 

total  X  hours  or  more  during  periods  t-l,t-2, . . .t-L  to  be  covered  in  period  t, 

it  is  clear  that  y.   will  be  correlated  with  h^  ,  for  t'  =  t-],t-2, . . .  t-L  and 

withy.  ,  for  t'  =  t-l,t-2, . . , ,t-L+l  and  for  t'  =  t+1,  t+2,,..,  t+L-1.  We 
It  

may  conveniently  display  the  pairwise  covariances  of  a  sequence  (y.  , 
T=l,2,ooc,t}  in  covariance  matrix  form: 


'11- 


^11 


IL 


^Ll 


o     « 


LL 


e     • 


L,2L-1 


Y  = 


*• 


\-L+l;t-2L+2'**   ''^t-L+l.t-L+l  *•'  "^t-L+l^t 


Y  Y 

t^t-L+1     ...   tt 


where  for  all  i,    Cov(y.  i,yi)=Y,   =q   ,-ppo  We  may  display 

XT  '    T  T       T  T      TT       T  T 

the  covariances  of  the  ys  and  hs  in  a  similar  fashion  as; 
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0       H^2      "l3   •  •  •  "l,L+l 


0       «23   •  •  •  "2,L+2 


L 
where  for  t  =  1,2, ...,t. 


"t-L,t 


"t-l,t 


.'\*  'N^ 


*N»  '^ 


-H  ,   =  E(e.  ,  y.  )  =  c  ,   for  7=r'+£,   1  <  i  <  L 

T  T        IT     it'      t  T  ~    - 


0  otherwise 
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2.4  Calculation  of  q 


T,T+i 


If  we  make  Assumption  IV: 


{e  ,    i=l,2^,.,,Nj  T=l,2, . . . ,t}is  a  double  sequence  of  mutually 

independent  identically  distributed  Normal  random  variables, 

2 

each  with  mean  0  and  variance  a   , 

then  for  each  i  and  for  1  <  i  <  L,  we  have  from  (15b)  and  (15d) 


Var(V,  )  =  iL-l)a 

XT  fe 


^"<W,^^^,)  =  Z  c^ 


Var  (r^^) 


£   a 


It  follows    that   the   (3  x  1)    random  vector 


~t 

i.    =  (V.      w.  r.  ) 

iT  IT  i,T+i  it' 

is  hfultinormal  with  mean  vector  (0  0  0)   and  covariance  matrix 


IL-i  0  0 

o.    0  i  0 

0  0  i 

Hence  for  all  i  the  (2  x  1)  random  vector 


*\^    *v. 


R  =  (R  R  ^ J  £  (r.  +  V.   V.   +  w,    ) 

—  T       -r+Ji  iT  IT  IT  i,T+i 


is  Multinormal  with  mean  vector  (0  0  )   and  covariance  matrix 


L-£ 


h~Z        L 


2 

1 

(1-|) 

e_ 

L 

_(l-t) 

1 

(22) 
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The  probability 


P(R     >  X   ,   R    .    >  X    , J    =  q        ^^ 

^     T  j'        T+i  T+i  T,T+i 


may  be  looked  up  in  tables  of  the  cumulative  unit-elliptical  bivariate 
Normal  function 


fN^^lJO,  |)  dz^  dz^ 


(23) 


h  k 


L(h,k,r)  E 

where  -1  <  r  <  1  and 
'l   r 
Z  = 

r   1 

tabulated  by  the  National  Bureau  of  Standards:  Tab les  of  the  Bivariate 
Normal  Distribution  Function  and  Related  Functions ,  Applied  Mathematics 
Series  50,  UoSoG.P.O.,  Washington  DoC,  1959, 

Notice  that  q      is  almost  in  proper  form  for  table  lookup, 

T,T+i 

We  need  only  observe  that 


where 


P(R  >  X  ,  R  .  >  X   )  =  P(z  >  z  ,  z  ^  >  z   ) 


T      t'  e  '   T     T  e  '   T+i      T+i  e  '  T+i      T+i  e  ' 


then  from  (23) 


P(R  >  X  ,  R    >  X   )  =  L( 

T      T     T+i      T+i 


nTlx   /lx 

T  T+i 


,1-1)     . 


However,  we  shall  use  the  computer  to  generate  values  of  q      needed  in  the 

^T,T+i 

course  of  the  simulation. 


-   31  - 
2.5     Calculation  of   c         ,. 

T  ,-T+i 


We  show  below  that 

if  we  make  Assumption  IV  stated  in  2.4;   For  i=l,2,.,„^N  and  T=l,2^..o,t 
the  e.   are  mutually  independent  identically  distributed  Normal  random 

IT  - 

2 
variables  with  mean  0  and  variance  o 

Proof;   From  (15a)  we  have 

00  00 


If 


and 


then 


Since 


=  /       I'  \ 


F^(z)  =F^(z|0,hp 


F^Cu)  =  F^(u|0,h^) 


*^T,T+^ 

(X     -U)  -00 


00 


/      ^  v^^v^- 


(^+i-> 


r      vf^^(v)  dv = + f^(/h^(x^^^-u)) 


we  have 


where 


Thus 
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=   e .  ^ b      /      e  Qu 

^/2n  v2jt        -00 


H  =  \\/i\-\)    and  A   =  X        h2/(h^+h^) 


T,-V+Z  — N*  T+i' 

N/2jt 
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3.   Outline  of  Simulation 

While  it  is  possible  to  describe  how  E(t]j^)  and  V(t]  )  behave 
as  a  function  of  (X,L)  quite  easily  using  (  10  ) ,  no  tractable  ana- 
lytical expression  for  the  probability  of  insolvency  given  n  within 
a  length  of  time  T  and  of  the  expected  waiting  time  to  insolvency 
given  Q,   appears  to  exist.   In  particular^  since  we  wish  to  describe 
how  this  probability  varies  with  changes  in  n,  monte  carlo  simula- 
tion is  a  convenient  and  flexible  method  to  use. 

The  steps  to  follow  in  simulating  values  of  P  (t)  and  of 
E  (t)  are  in  a  broad  outline  listed  below.   For  each  value  of  f^ : 


(1)  Calculate  q    ,.,  P  >  c    ,,  for  1  <  i  <  L 

TfT+l'      t'  T,r+£  -       - 

and  1  <  T ,7+Z  <   2t  as  described  in  2.4  and  2,5, 
—     '  —   o 

(2)  Use  q    ,  . ,   P    ,    c.         ..*  k»  and  r  to  calcu- 
late  elements  of  Z  and  of  6  as  shown  in 

(  6  ), 

(3)  Calculate  S  as  described  in  3.2. 

(4)  Carry  out  the  simulation  routine  flow  dia- 
grammed in  section  3,2. 

(5)  Estimate  P  (t)  and  E  (t)  and  calculate  the 
variances  of  these  estimates. 


The  simulation  routine  consists  of  repeating  the  following 
steps  a  large  number  of  times,  say  R  times: 

(i)   Given  t,  generate  a  sequence  of  t  standard- 
ized Normal  random  deviates,  {u  ,t=1,2, , , , ,t Jj 
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(ii)   use  6   from  (  6  )  and  a  matrix  S  such  that 

-.t"*^  ^        /     ~    .t  . 

S   ^    -  L   to   transform,  u^s  (u,^u^.o»u^)   into 


-t-  ^1^  2" 
to  compute  S  in  3  2  i „ 


a  simulated  observation  of  6  »   (We  show  how 


(iii)   check  to  see  if  there  is  a  iv  £  ^^  1  <  t  <  t, 
If  so^  record  that  ruin  occurred  on  this 
particular  replication.   Record  also  the 
smallest  t  for  which  6  <  0  if  ruin  did 


occur 


(iv)   Repeat  (i);,  (ii)>  and  (iii)  R  times, 


To  estimate  P  (t)  from,  the  results  of  the  simulation  routine 


we  let 

1         insolvency  occurs  on  the  ith  replication 

X.  =J     if 

^   1 

otherwise 

and  regard  (x.  ^i:rl  ^2^,  o » ,  ,R}  as  a  sequence  of  (independent)  observa- 
tions generated  by  a  Berno'Jii.lli  process  with  parameter  P  (t)  <.   An 
unbiased  estimate  of  P  (t)  is 

-   1   ^ 

R  1.-1   1 

and  an  unbiased  estimate  of  the  variance  V(x)  cf  TT  is 


~     1   ^         2 

v^^^  =^  i::r  ih  ^\'  ^) 


* 

We  will  use  a  sequence  of  such  t's  to  estimate  the  expected 
W'_itin,^  time  F._^i't)  to  insolvency  conditional  on  ruin  occurring 
for  some  1  <  -■-'-<'  r  . 
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We  estimate  E  (t)  and  the  variance  of  this  estimate  in  a 

n 

similar  fashion;   let 

insolvency  occur-:  at  I  <  •^.  <  t 


w.  =  '^      if 


otherwise 


and  regard  (w.  ^i=l ^,2;  . ,  o^R}  as  a  sequence  of  independent^  identi- 
cally distributed  random  variables.   Then  an  unbiased  estimate  of 


E  (t)  is 


w^J^w.  (26) 


and  an  unbiased  estimate  of  the  variance  V(v^)  of  w  is 

~     1    R        9 
V(w)  =  -^  .Ej^  (w.-w)  ,  (27  ) 
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3  c  2   Calculation  of  S 

In  order  to  carry  out  the  simulation,  we  nwst  find  S   from  E^  which  in 

turn  is  given  in  terms  of  p^  r,  k^  H  and  y  in  (21b)  „   Since  the  order   t  of 

E  required  during  the  course  of  the  simulation  is  in  general  very  large 

(e„g„  t  =  500  is  not  unreasonable)  we  will  show  how  cyclicality  of  elements  of 

E  may  be  exploited  to  make  calculation  of  S   not  only  computationally 

feasible^  but  also  analytically  simple,, 
the  application  of 
In  Part  II  of  this  paper  we  show  that  the  covariance  terms  a      ,      of 

TT 

E  are  periodic  in  the  sense  that  there  is  an  integer  t   such  that 
=  o 


TT       T+KT  ,T     +KT 
O         O 


(  28) 


for  all   K=0,l,2, 


sue 


h  that   1  <  T+KT  <  t   and   1  <  t'+kt  <  t.  When 


o  — 


o  — 


the  order   t  of  E   is  a  multiple  of  t   (28)  allows  us  to  write  E  in 

=  o  = 

the  form 


E  = 


(29) 


C    B    C 


C    B 
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where 


T   J.   0     c 

o 


It 


and 


C  - 


'l,T  +1 

'   o 


't  ,r  +1 
o     o 


1,2t 


(30) 


2t  ,2t  * 

O     O  I 


We  show  below  that  S  of  any  order  that  is  a  multiple  of  t  has  these 

=  o 

properties  provided  (  28  )  is  true; 

(i)   To  find  S  we  must  diagonalize  at  most  four  positive 

definite  symmetric  matrices--two  of  order  2t  and  one 

of  order  j    «   (When  the  order  of  S  is  an  evefi 

multiple  of  t  ,    only  two  matrices~of  order  2t  and  one 

of  order  t  need  to  be  diagonalized,) 
o 

t  i  i 

(ii)   The  matrix  S  =  (U  M)  a^  ,  where  UjM  and  A^  are  defined 

in  terms  of~matri;ces  derived  in  the~course  of  the 

diagonalizations  mentioned  in  (i). 

Exhibit  2  outlines  an  easy  method  for  calculating  all  but  the  last  2t 

o 

rows  of  S   in  terms  of  B.>  C^  the  orthogonal  matrices  P  and  Q  that 
diagonalize  B  and  ^   respectively,  and  the  diagonal  matrices  A  arid  D 
defined  in  (3i)    and  (  38  )  below. 

PROOF g  We  shall  deal  with  the  case  when  the  order  of  Z  is  an  odd  multiple 
of  T  o   The  modification  necessary  when  the  order  is  an  even  multiple  of 

T  will  become  evident „  Si. ice  the  matrix 
o 


E  s 


t 
C    ii 


(31) 
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is  positive  definite  symmetric^  there  is  an  orthogonal  transformation  Q  of 

dimension  (2^     x  2t  )  which  diagonalizes  E; 
o     o  °  = 


Q  E  Q 


^2- 


-  A 


(32) 


Partitioning  A  into  (t  x  t  )  matrices  we  may  write 
=        o    o 


A  = 


en 


§22 


(33) 


Defining  the  (t  x  t)  orthogonal  matrix  M  composed  of  Q's  plus  a  (t  x  t  ) 
identity  matrix  down  the  diagonal  and  zeros  elsewhere. 


M  = 


(34) 


it  follows  that 


Ta,, 

I  =11 


^22 


I 


M  Z  M  =: 


t 

L 


ei] 


^2 

t 


11 


^22    e 


(35) 


A3 


Letting 


V  = 


^22 


^11 


^J 


and   V 


^22   £ 


L= 


(36) 


we  may  write 

A 

M  L  M"  == 


11 


V  C 

~   V 


V 
=o_j 


(37) 


Since  V  is  also  positive  definite  synmietric  there  is  a  (2t  x  2t  )  orthogonal 
=  ■   o     o 

matrix  P  which  diagonalizes  Vj    that  is. 


P  V  P   = 


=  D   . 


(38) 


19 


2t 


Similarly,  let  P  be  the  orthogonal  matrix  which  diagonalizes  V  ,    and  D   the 
corresponding  diagonal  matrix. 

Defining 


U  = 


(39) 


P 
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we  may  write 


U  M  E  (U  M)   =  A 


where  A  is  the  diagonal  matrix  whose  diagonal  elements  are  the  characteristic 

roots  of  Z. 

Thus 


Z  =  [(U  M)''  A^^l  [(U  M)''  A^]*" 


(40) 


or 

E  =  S  S* 
where 


S  =  (U  M)  A= 


(41) 


We  may  use  (  34  )  and  (  39  )  to  determine  S  in  terms  of  partitioned 
elements  of  Q  and  P.   Partition  both  P  and  0  into  elements  composed  of 


(t  X  t  )  matrices: 
o    o 


P  = 


P     P 
=11   =121 


=21 


=22j 


p(°) 

=11 

Jo) 
=12 

P    = 

=o 

p(o) 
_=21 

p(°) 
=22_ 

Sii   ei2 


Q  = 


iS21   522 


(42) 


Then 
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U  M  = 


111     In 


111     111 


(o) 

p 
=11 


p(o) 
=21 


p(o) 
=12 


,(°) 
=22_ 


2ll        111 


111       222 


Sll        9l2 


?21        222 


^11 
<ill521> 


2l2 

(PnQ22) 


(|2l52l)  (l21?22) 


(luiiP 

(|22eil) 


<?ll22l) 


<l2ie2l) 


(Il25l2) 


(P222l2> 


(Ell222) 


(?2l222> 


(!l22ll>         ^!l29l2) 


(|222ll>         (|222l2) 


(C) 


<?2l22l)  (l2l222^        I22 


(43) 


Partitioning  D   and  D      into   elements   composed  of    (t     x  t   )    matrices   we  have 
=  =0  ^  00 


Ell      2 


2  ^22 


D     = 

=0 


(o) 

?u       2 


0  Di°> 

=22 


(44) 


Remembering  that 
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r  A. 


A     = 


11 


So 


(45) 


We  obtain  from  (    41  ), 


s'  = 


Z        0 
G       W 

G     W 


G     W 


G     W 
=o   =o 


where 


I  =   (Ml§ll  ^Il9l2) 


D      = 

siO 


=  11 


111, 


G   = 


flllulll^  (i!l?llS22^ 

>^El2l2l52l)  ^?t2E2lS22) 


G 

=0 


^Eii^lu  92,^ 

.  h        (o)    (o) 

^=22  =21  =21'' 
u.       o 


^=11   =11   522'' 
o 

i        (o) 

(522   ?21  222^ 


W  = 


(e!lll2eil>  (H!i^122i2) 

\:flllllll{>  ^Et2l222l2> 


w 

=o 


D*     P(°^ 

=11   =12 

o 

=12  =22  J 


-  47 


This  completes  the  proof  of  properties  (i)  and  (ii)  when  the  order  of 

S  is  an  odd  multiple  of  t  .  When  the  order  of  S  is  an  even  multiple  of 
=  o  - 

T  V  =  A   which  is  already  a  diagonal  matrix.   Hence  the  second  part  of 
o  =o   =22 

* 
(i)  must  be  true. 


In  the  course  of  the  simulation  described  in  Part  II,  we  will  assume 
that  the  order  of  E  and  hence  of  S  is  so  large  as  to  obviate  the  need  for 
diagonalizing  V  even  if  the  order  of  S  is  odd.   Hence  we  have  omitted 
description  of  it  in  Exhibit  2, 
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EXHIBIT    2 
CALCULATION  OF  S 


Find  an  orthogonal  matrix  Q  e 


Sii     9i2 

921   922. 


of  order  2t   that  diagonalizas 
o 


E  = 


B   C 


t 
C   B 


That  is^  a  Q  such  that 


Sii    9i2 


921   922 


B    C 


C^   B 


r 

9u   921 


t      t 

12  522 


^11   9 


2     ^22 


where  A   and  A   are  diagonal  matrices  of  order  t  . 


Find  an  orthogonal  matrix  P  = 


F     P 
=11   =12 


P     P 
=21   =22 


that  diagonalizes 


V  = 


^22   9 


.9        en 


This  is,  a  P  such  that 


r 


p    P 

=11   =12 


P     P 
=21   =22 


^: 


A, 


11 


t    t 
P    P 
=11   =21 


t    t 
P    P 

=12  =22 


Ell     9 


9     922 


where  D   and  D„„  are  diagonal  matrices  of  order  t  < 
=11     =22  o 
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(EXHIBIT    2  CONTINUED) 


Calculate 

the  matrices 

Z   = 

<e!iEii  4si2> 

klvlirhi^ 

^gi^.irS22>~ 

G  = 

(?22"l2rS21> 

(Et2'?2rS22> 

^Eu°El2»Sll) 

(iir£i2°8i2>' 

W  = 

sii'hrhi^ 

^4"E22°Sl2>. 

> 

f 

Write  all 

but   the   last   2t     rows 
o 

of  S^   as 

Z 

0 

G 

W 

l'- 

G     W 
G     W 

o            o 

o             o 

- 

O                0 

